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Abstract 

We consider right-handed neutrinos propagating in 5 (large) extra dimensions, whose only cou- 
pling to Standard Model fields is the Yukawa coupling to the left-handed neutrino and the Higgs 
boson. These theories are attractive as they can explain the smallness of the neutrino mass, as 
has already been shown. We show that if 5 is bigger than two, there are strong constraints on the 
radius of the extra dimensions, resulting from the experimental limit on the probability of an active 
state to mix into the large number of sterile Kaluza-Klein states of the bulk neutrino. We also 
calculate the bounds on the radius resulting from requiring that perturbative unitarity be valid in 
the theory, in an imagined Higgs-Higgs scattering channel. 
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I. INTRODUCTION 



The Standard Model (SM) of high energy physics suffers from the gauge hierarchy prob- 
lem, which is the fine tuning required to maintain a low electroweak scale (Mew ~ 10 3 GeV) 
in the presence of another seemingly fundamental scale, the Planck scale (the scale of gravity, 
M p i ~ 10 19 GeV). Supersymmetry, technicolor and more recently extra (space) dimensions 
have been proposed to address the hierarchy problem. 

Recent neutrino oscillation experiments have suggested a non-zero neutrino mass, with 
the best fit values of the mass differences and mixing angles given by 

AmLar = 7 x 1(T 5 eV 2 , tan 2 9 Solar = 0.4 , 

Am 2 Atm = 2.5 x KT 3 eV 2 , tan 2 9 Atm = 1 . (1) 

The oscillation between the three active flavors, v e ^v^v T , in the SM, accommodates this 
satisfactorily, with Amg olar = \m\ — m|| and Am\ tra = (m 2 , — m||, the rrii being the physical 
masses. If the m ; are also assumed to be of the same order of magnitude as the mass 
differences, it is quite challenging to explain why it is that the neutrinos are so light compared 
to the other leptons. 

It has been shown 00 that if there ate other Large Extra Dimensions (LED) in adchtron 
to our usual four space-time dimensions, we could potentially solve the gauge hierarchy 
problem. It was then pointed out 0, 0, that the smallness of the neutrino mass is 
naturally explained 2 if right-handed neutrinos that propagate in some 5 number of these 
extra dimensions are introduced. We will refer to such neutrinos, which are SM gauge 
singlets, as "bulk neutrinos" , as is the usual practice. Various aspects of theories with bulk 
neutrinos have been analyzed in Ref. ^jj. 

If bulk right-handed neutrinos are responsible for the smallness of the neutrino mass, we 
ask in this paper, what the constraints on such a model might be. We will find that the 
experimental constraints on the probability of an active neutrino to oscillate into sterile bulk 
neutrinos, can give us a lower bound on the inverse radius (1/-R) of such extra dimensions, 
especially if 5 > 2. We will also require that perturbative unitarity in two-to-two scattering 
of Higgs bosons be preserved in the theory. We will include the tower of Kaluza-Klein states 
as intermediate states, and again find that for 5 > 2 it results in a strong bound on 1/R. 



We will take the view, as in Ref. 



12|, that the standard three-active-flavor oscillation 
explains the data in Eq. and that the mixing to sterile bulk neutrinos are small enough 
to evade experimental constraints. Here we extend the analysis to 5 > 1 and show that 



1 In this work, we will not address the LSND result 

2 We note here that the conventional see-saw mechanism to explain the smallness of the neutrino mass is 
equally appealing, but we will not consider it in this work. 
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strong bounds on 1/R can result in that case. However, a precise calculation of the bound 
will not be possible when 8 > 2 due to sensitivity on the cutoff scale, implying a dependence 
on how one completes the extra dimensional (effective) theory that we will work with. As 



pointed out in Ref. |l2j an alternative approach [8|, [lOj], wherein Eq. (JTJ) is explained by 
the oscillation of the active species predominantly into sterile bulk neutrinos, appears to be 
disfavored by the Sudbury Neutrino Observatory (SNO) [2| neutral current data. 

In order to constrain theories with neutrinos in the bulk, other processes have also been 
considered in the literature 0] including muon lifetime, pion decay, flavor violation, beta 
decay in nuclei, muon g — 2, supernova energy loss and big bang nucleosynthesis. Some 
phenomenological implications are also considered in Ref. 14j . 

The rest of the paper is organized as follows. We will introduce the extra dimensional 
theory with bulk neutrinos in Sec. |Hl write down the equivalent four dimensional Kaluza- 
Klein theory with particular focus on the interaction of the right-handed bulk neutrino with 
the Higgs field and the left-handed neutrino, and diagonalize the neutrino mass matrix. In 
Sec. IIII1 we will find the bounds on the radius of the extra dimensions from limits on neutrino 
oscillation into sterile states, and in Sec. IIVI from unitarity considerations in Higgs-Higgs 
scattering. Our conclusions are given in Sec. We will present an alternate approach to 
the diagonalization of the neutrino mass matrix in Appendix El and give in Appendix [B] 
some detailed formulas that we use in deriving the unitarity bound. 



II. RIGHT-HANDED NEUTRINOS IN EXTRA DIMENSIONS 

To address the gauge hierarchy problem, Arkani-Hamed, Dimopoulos and Dvali (ADD) [f| 
postulate that the Standard Model (SM) fields are confined to a 4 dimensional (4-D) sub- 
space (brane) in an extra dimensional world of 4 + n total dimensions. ADD take the view 
that the only fundamental scale in nature is M*, which is of the order of Mew-, an d the 
apparent 4-D gravity scale (M p i) is then given by 

Ml = M* +n V n , (2) 

where, V n is the volume of the (compact) extra dimensional space. In the simple case of 
each of the compact extra dimensions being of equal radius R', we have V n ~ R' n . Thus 
ADD argue that M v \ appears to be a large scale from a 4 dimensional perspective simply 
because the volume V n is large. In other words, the explanation of why M p i is large is recast 
to stabilizing R' at a large value, so that V n is large. 

It should be pointed out that for a given M*, if the n compact dimensions have equal 
radii R', Eq. (J2j) implies a particular value of R'. However, if it happens that there are two 
sets of compact extra dimensions of unequal size, 8 of them (8 < n) with radius R, and the 
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other (n — 5) with radius R', then we have in this case, V n ~ R' {n - 5) R 5 . We can in this case 
think of R as an independent variable with R 1 being determined by Eq. (j2J. 

We consider the ADD framework, to which is added three (one for each generation) bulk 
fermions, ^ a (x^,y), that propagate in 4+5 dimensions (5 of them compact with radius R), 
where the indices a, (3 — (1, 2, 3) denote the three generations, and y stands for {y 1 , ...,y s }. 
Since the higher dimensional Lorentz invariance is reduced to our 4-dimensional invariance 
by the presence of the brane (where SM fields are localized), we are interested only in 
keeping track of the 4-dimensional Lorentz transformation property of \I/ along the direction 
of the brane. We consider the situation where the SM fields couple only to 2 components 
of denoted as ip^ix^^y), and transforming as a right-handed 2-component Weyl spinor 
under 4-D Lorentz transformations. For example, for 5 = 1, we denote a 5 dimensional (4 
component) spinor as 

where the L and R subscripts make explicit the four dimensional Lorentz property, and in 
particular, ip® R each transforms as a (2 component) Weyl spinor. In the following, we will 
keep the formalism general enough to include arbitrary number of extra dimensions. 
We can split the Lagrangian into a bulk piece and a brane piece, 

S = J d 4 X d 5 y [£ B ulk + 5 (y) £ Brane] • (4) 

^Buik contains the Einstein-Hilbert bulk gravity term (which we will not show explicitly, but 
can be found, for example, in Ref. |15j). the kinetic energy term for the bulk neutrino field 
^/(x^ 1 , y) and in general, a bulk Majorana mass term for which for simplicity we will omit 
(see Ref. |8| for implications of a nonzero bulk Majorana mass). Z^Brane contains the SM 
Lagrangian plus an interaction term between SM fields and ipR, 

£ Bu ik d * a ir M D M * a , 

£ B rane D ^SM ~ ( — = h tjj^ + h.C. ) , (5) 



where, is an 0(1) Yukawa coupling constant. It should be kept in mind that ipR is a 
function of (x M , y) whereas the SM fields are functions of x M only. The index M runs over 
{x»,y}. 

We can perform a Kaluza-Klein (KK) expansion of the 4+5 dimensional theory and obtain 



4 



an equivalent 4-D theory by writing, 



^sk^E^)^. ( 6 ) 

n 

where, n = (ni, ■■■,ng) is a vector in "number space", ip^ are the KK modes and fn(y) is a 
complete set over y. A similar expansion is made for i/j^. To reduce clutter, we will simply 
write n and y for n and y respectively. We will use the notation n = (0, 1, ...) and n = (1, ...) 
(n excludes 0). /„ is an orthonormal set, 

r R d 5 y r n (y)f m (y) = 5 nm , (7) 
Jo 

and a convenient choice is, 

Uv) = = % (8) 

g ~ 1 R s v 5 



with Sa-i the surface "area" of a unit sphere in § dimensions, and V$ = ^f^-R 6 is the volume 
of the extra dimensional space. 
We define the fields 3 |l2f, 

s _L(^) + ^ ( - fi) ) . 

We substitute the KK expansion for the bulk fields ipR an d V'x m t° Eq. © to get the 
equivalent 4-D theory, 



a=l h 



a,/3=l L 



^ + ^72^^ + /i.e. 



(9) 

where, \n\ = v/nf + ... + njj. We note here that in there is a tower of KK states 
(y^ i v ^) with Dirac masses approximately equal to \n\/R. Henceforth, we will assume 
that unless noted otherwise, repeated generation indices a,j3,i,j, and KK indices n,n,m, 
are summed over. 



The oilier linear comlimat ions (i>^"^ ~ ^Pr and y^"^ ~ ^i} are decoupled from ilnSM li< ■!< Is 



and we will not consider them further. Also, with orbifold compactification, we can project out ipi,^ so 
that it is excluded from the particle spectrum. 
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TABLE I: Estimate of A^. (M* = 1 TeV) 



<5 






1 


1 




2 


1 


0.1 


3 


10 3 


1 



With 577(2) broken by the Higgs mechanism, by the Higgs field acquiring a vacuum 
expectation value (VEV), (h) = v, we have the neutrino masses given by, 

, (10) 

and m„ can be much smaller than M^w ~ t>, if K5 is large. We show in Tableman estimate 
of the order of magnitude of the A^ that one needs in order to get neutrino masses of 
O(10~~ 2 ) eV. The 1/R values are chosen anticipating the constraint from neutrino oscillation 
that we will derive in Sec. IHII We note that 5=1, though the case most studied, requires 
an unnaturally small A. We will therefore include 5 > 1 in our analysis. 



A. Truncation of the KK sum 

Though hi can in principle go up to 00, leading to an infinite tower of Dirac states with 
masses \h\/R, we take the view that this extra dimensional field theory description is valid 
up to the cutoff scale M*, and therefore truncate the hi such that the highest KK mass is 
M*. We define N to be N/R = M*. For 5 > 1, the state with mass ^ can be degenerate, 
and we denote the degeneracy at the h th level by d^. (Strictly speaking, we should denote 
this as d\n\, but we will just write this as d^.) For example, for 5 = 3, the state with mass 
1/R has d\ = 3, corresponding to {hi,h 2 ,hz) — > (1,0,0), (0,1,0) and (0,0,1), all of which 
have the same mass. For large \h\, the leading power dependence of da in 5 extra dimensions 
is given by d^ = Cf^h^^ 1 , where the c„ are 0(1) numbers. We define do = 1. 

For large \h\, we can think of the hi as a continuum and the leading behavior is given by 
the surface of the (5 — l)-sphere of radius \h\ in number space, 

dh ~ *S' ( 5_i |ri|^ — 1 (in 5 dimensions). (11) 
4 =1 = 1, 4= 2 ~ 27m, 4= 3 ~ 47rn 2 . 

For example, for 5 = 3, d^ ~ 47r|n| 2 , which is the surface of the 2-sphere with radius |n|. 



6 



Thus, N is the radius of the biggest sphere in {rii} space such that N/R = M*. We will 
often use the continuum approximation for estimating various quantities. 

The sum over the KK states of certain quantities can be divergent and can depend on N. 
We will elaborate more on this later, and we will see that the probability of active neutrinos 
oscillating into heavier sterile KK states can depend on N, especially strongly for 5 > 3. 



B. Mass Matrix Diagonalization 



At the outset, we can make the Yukawa coupling in Eq. © diagonal in generation space 
(a, (3) with the rotations [l^j ]. 



V 



R 



r v 



L J 
li(h) 
R J 



"R = (OX , 



r ai >i{n) 



V 



a(ii) 



R 



{r™)*e 11 



L 

H 

R i 



where the unitary matrices I and r are chosen to diagonalize m^, so that, (r m )* m^f (i™) = 
m\,8 %3 . Similarly, l e and r e are chosen to diagonalize the electron- type mass matrix. After 
these rotations, Eq. Q becomes, 



r (4) r \ n \ ( H(n) H(h) . , 

C ( > = C SM --^\ U R V L + h - c - 



;i2) 



In the usual way, the charged current interactions now become proportional to the MNS 
matrix 

V M ns = 1\1. (13) 



By defining the Dirac spinors, 



v 



v 



L 

4. 



V'^ EE 



V 



V 



«(1) 

L 

H{1) 
R 



. V 



H(n) — 



V 



V 



L 

li(n) 
R 



(14) 



we can rewrite Eq. (jT2j) as, 



£(4) = r _ M^ftJ^W _ IhL 



R 



hu H u H + V2h v'^Plv' 1 + h.c. 



(15) 



where P L{R) = (1=F7b)- 

To reduce clutter, in the rest of this section, we will suppress writing the generation index 
i and the prime on the fields. Thus for each i = 1,2,3, the mass matrix in KK space can be 
diagonalized independently using the procedure described next. 
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For each i, the neutrino mass term obtained by setting (h) = v is, 

with (^d) T = iy v^> ■■■ ...) and the mass matrix is given by, 

/ m v \f2m v P R ■ ■ ■ V2m u P R . . 
V2m u P L 



(16) 



M 



D 



■ \ 



1_ 

R 



V R Ja^xclf, 



V • 

We keep in mind that the n th state is degenerate with degeneracy d fl . 
We define 

i = m u R, ^ = -f-. 

n 



(17) 



\M i 

R / 



(18) 



We will restrict ourselves to the situation when J2h^fi^n <C 1 (for all S), since, as we will 
show later, this is the condition implied by experimental data on oscillation into sterile 

states. 

We can diagonalize M.d by two unitary matrices L and R so that 



v D = {LP L + RP R )v D , 



(19) 



where v D is the mass eigenvector. We perform the diagonalization perturbatively, and to 
0(£ 2 ), the L and R are, 

/(1-E^n) V^l ... >/2& • • V^N \ 

-vki (i-£i 2 ) 



(1 £l)dftXd« 



(20) 



A? 



/ 1 ■ ■ 

-V2g 1 



R 



V2ti ■ ■ V2&\ 



(1) 



dfi, x d n 



(21) 



The mass eigenvalues are the diagonal elements of R)M.dL. (See Appendix El for an 
alternate approach to diagonalization and for details on obtaining the mass eigenvalues.) 
We find the lowest mass eigenvalue, 



rrv ' m m,. 



\n\ 



(22) 



Given that we take £ to be small enough so that J^£|<in <S 1, we have « m v . We 
note that henceforth, when we write m, we mean m^°\ The heavier mass eigenvalues 
at the n th level, with gL states, are 



rn 



in) 




[l + (£ 2 + £ 4 )dft] for 1 state, 

for the other (dg — 1) states. 



(23) 



In summary, reintroducing the generation index i and the KK index n — ► {0,n}; n 
1, iV, we can write the flavor state in terms of the mass eigenstates as 



,a _ jai t On 



(24) 



We can write a similar expression for the right-handed neutrino, but we will not need that. 



III. NEUTRINO OSCILLATION CONSTRAINTS 

Neutrino oscillation experiments measure the probability of producing an active neutrino 
flavor a, and detecting an active flavor (3 a distance L away. Recent experiments, particularly 
Super-Kamiokande and SNO, indicate that the oscillation among the three active species 
v e , v T is the best fit to data. We will denote the probability of an active species v a 
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oscillating into another active species vp after a distance L as P Va ^, V(j . The mixing of an 
active species into a sterile species u s , P Ua ^ Vs = 1 — J2p Pv a -*v -, is furthermore, strongly 
constrained by CHOOZ and by fits to the atmospheric oscillation data. We will thus 
assume that the standard three active species oscillation with the usual MNS matrix provides 
the best fit to data, and the oscillation into the heavier KK states (sterile states), would 
have to be small enough not to violate the constraints on P Ua _> Us . This was the view taken 
also in Ref. jl^ . where it is noted that the constraints on the mixing into sterile states are 
the following (at the 90% C.L.): 



CHOOZ 
Atmospheric v 



P v ^ Vb < 0.058, 



< 0.17 



\[P v ^ v , + P Ve ^ Vm ] <0.39 



(25) 
(26) 



We have chosen these bounds as they place the strongest constraints on our model. 

The probability for an active species known to be v a at t — 0, to oscillate into the active 
species vp after time t (after traveling a length L), is determined by the Hamiltonian H, and 
is given by the time evolution operator e~ tHt . (For a concise account on neutrino oscillation, 
see for example Ref. 0]) Therefore, 

2 



P 



Using Eq. (j23| to write z/£ in terms of the energy eigenstates u^ n \ and since 

H z>i (n) = z>i (n) , 



(27) 



(28) 



we get, 



l^r\L^\ 2 d n e 



(29) 



where Lf 1 is the first row in Eq. (|20|) and d Q = 1. For a neutrino beam with energy E v and 
momentum p u , in the relativistic limit, we have, E 1 ^ « p u + m\ /(2E V ), and therefore, 



P 



(30) 



m 



( n \ the mass eigenvalues, are as given in Eqs. (|22|) and ()23j). From Eq. (j2*U)l we have 



- oo 



z — ' rr 



L 



On 



V2 



n 



(31) 



which implies to 0(£ 2 



/3 



AE U R 2 J ' 



(32) 
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where we assume that | > mj, which we will see soon is required to satisfy the experimental 
constraints. Hence, the probability for an active state to oscillate into sterile states is 



Before we can use Eq. to compute the bounds on 1/R, we have to specify the 
rrii and the l ai . We take the rrii given by Eq. (|22|) to be consistent with the mass dif- 
ferences extracted from standard fits to the solar and atmospheric neutrino oscillation 
data given in Eq. So far we have precise information only on the mass differences 

from the oscillation data, with upper limits on the masses otherwise. The limits are 
(m Ue , m v , m„ T ) < (3 eV, 0.19 MeV, 18.2 MeV) from final state lepton spectrum, and, 

the sum of the neutrino masses . m v . < few eV, from cosmology (for reviews, see Ref. j^J). 
It is usual to consider the following three mass schemes, all of which are consistent with 



(i) Normal Hierarchy: m\ ~ 0, m 2 pa 0.008 eV and m 3 pa 0.05 eV. 

(ii) Inverted Hierarchy: m\, rri2 ~ 0.05 eV and pa 0. 

(iii) Degenerate: m\, 1712,1713 all at some mass scale less than the limits given above. Here, 
to illustrate the character of the bounds, we take the three masses to be around 1 eV, 
an arbitrary choice. 

We can work in the basis in which the charged lepton mass matrix is diagonal. In this case 
we have from Eq. (|13|h Vmns — I- The bounds we derive are not very sensitive to the precise 
values chosen for I. Based on the Vmns obtained from a global fit to the oscillation data, 
we take the I to be such that the solar neutrino mixing angle between u e & is given by 
tan 2 9 efl = 0.4, and the atmospheric oscillation mixing between & v T is maximal. This 
implies 



Given a particular mass scheme, we can then derive the bounds on 1/R using Eq. (|3*3*|) . 
In this equation we note that for large KK masses, the argument of the sine function is 
large, leading to a rapid oscillation with distance. We therefore use the average value 




(33) 



Eq. ©: 




(34) 



sin 2 (4^2) ) = 1/2, and require that the experimental constraint be satisfied by 





(35) 
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TABLE II: Lower bound on -s (eV) for right-handed neutrinos in 5 — 1,2,3 extra dimensions. 





Normal 


Inverted 


Degenerate (m « 1 eV) 


CHOOZ 


Atm 


CHOOZ 


Atm 


CHOOZ 


Atm 


5 = 1 


0.03 


0.15 


0.5 


0.13 


10.6 


4.1 


8 = 2 


0.32 


1.5 


5.3 


1.3 


100 


41.7 


5 = 3 


2.4 x 10 3 


5.6 x 10 3 


1.2 x 10 4 


4.9 x 10 3 


10 5 


5 x 10 4 



We remind the reader that in the continuum approximation, the da is given in Eq. (jlljl . The 
bounds on 1/R obtained by requiring Eq. (}35|) to satisfy the limits in Eqs. (}25|) and (J26|) are 
summarized in Table |H] For the lighter KK states, with 1/R ~ 1 eV, as could be the case 
for 5 = 1,2, the argument of the sine function in Eq. ()33|) may not be large and a numerical 
calculation may be necessary. However, since the number of such states are small, we expect 
that the bounds we found using the average value is a good approximation. Ref. has 
presented the bounds for 5 = 1 by performing an exact numerical calculation, and we find 
that our bounds are quite similar to theirs. 

We would like to stress here that for 6 > 1, the sum over the KK states, for example 
in Eq. (J33j) . is not convergent, given that the degeneracy at the n th level is d n ~ n 5 ^ 1 . To 
illustrate this further, considering 5 = 3 with d n ~ Aim 2 , we get equal probability, for all n, 
for the active neutrino to mix to the n th KK level, after including the degeneracy. Therefore 
our estimates are valid only up to 0(1) coefficients since we are sensitive to the manner 
in which our four dimensional KK description is completed to a more fundamental theory 
at the scale M*, the cutoff of our theory. 4 An additional source of uncertainty comes from 
using the continuum approximation rather than performing the discrete sum over n. 

One might wonder if in a certain high energy completion (or if one includes higher di- 
mension operators in the effective theory), it is possible that the off-diagonal entries in the 
neutrino mass matrix, Eq. (fTTjl . could be vanishingly small compared to the 00 entry, so 
that the above mentioned oscillation constraints could be severely weakened. Since it is the 
same operator that generates the 00 and the off-diagonal entries, c.f. Eqs. (0) and this 
would mean that the 00 term also has to be vanishingly small, which would conflict with 
the experimental constraint in Eq. (PJ. We thus conclude that the oscillation bounds on 
1/R that we estimated above are general to the class of theories, in which only right-handed 
neutrinos can propagate in bulk, though they depend to some level on how the theory is 



For a related discussion in the graviton sector, see Ref. [21| 
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completed at the scale M*. For 5 > 3 the active state would oscillate mostly to the heaviest 
states and we would not be able to reliably estimate the oscillation probability using the 
effective theory that we are working with. We thus restrict out conclusions to 5 < 3. 



IV. UNITARITY CONSTRAINTS 

We would like our extra dimensional theory to be unitary at each order in perturbation 
in order to have the ability to make reliable predictions. The amplitude for certain processes 
that receive contributions from a tower of KK intermediate states can grow with the center- 
of-mass (cm.) collision energy yfs. For such processes, the requirement of perturbative 
unitarity can be particularly stringent as yfs increases. 

For instance, the number of KK states accessible at a collision energy yfs is of the order 
of (y/sR) s , which can be large. For the amplitude of a scattering process that grows with 
the number of KK states in the intermediate state, requiring that unitarity be maintained 



for ^/s < M* can lead to a lower bound on 1/R. As pointed out in Sec. Ill Al the sum over 
the KK states can be divergent, and therefore, we cut the sum off at N s = ^JsR. 

We consider the imagined scattering of a pair of Higgs bosons hh — ► hh, as shown in Fig.[IJ 
to study the unitarity constraints on 1/R as a function of yfs. The two incoming (outgoing) 
momenta of the Higgs bosons are labeled by p\ and p 2 (ps and p^j. With p\ +p 2 = P3+P4, we 
define the Mandelstam variables in the usual way, s = 2pi -p2 + 2m|, t = —1p\ •£>4 + 2m^ and 
u = —s — t + Aml = —2pi -£> 3 + 2m^. The kinematics of 2 — > 2 scattering are simply described 
by s and one scattering angle, 9, in terms of which t = — (l — Am 2 h / s) (1 — cos9) and 

u = — — (l — Am\/ s) (l + cos6>). The amplitude A can be expanded in Legendre polynomials 

as 

00 

A(hh -> hh) = (32vr) ^(21 + l)Pi(cos9) a, . (36) 

1=0 

Given the scattering amplitude A, the partial wave amplitude ai are determined by 

• 1 



ai - 

QAtt 



^- J dcos# Pi (cos 9) A(hh -> hh) . (37) 

For the process we are considering, for unitarity to be satisfied, it is sufficient that 

Im(a ) < 1 , (38) 

where a is the S-wave amplitude. We will thus focus on calculating Im(^4) in order to apply 
this bound. 
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FIG. 1: Diagrams that dominantly contribute to Im(^4). 



The dominant contribution to the imaginary part of the scattering amplitude A can be 
written as 



Im (A) = lmlj2 M nm d n d m j 

\ nm / 



(39) 



where the Feynman graphs for A4 nm are shown in Fig. in which we have shown only the 
diagrams that contribute dominantly to Im(^4). 

Before detailing the full calculation, we present an order of magnitude estimate of Im(^4). 
Due to the unitarity of the S-matrix, Im(^4) is given by the sum of the cut diagrams. We 
obtain an estimate by ignoring the KK masses in the intermediate state, and furthermore, 
by using the continuum approximation, c.f. Eq. (jllj) . We define N s so that the highest KK 
state accessible to the cut diagram is given by ^ = y/s, and we get 

i N ' 4 

Im <-4) ~ ^£(v) (ft-iInT) (ft-M*- 1 ) , 



1 fm v \^ 1 ( S S -i 
167T \ v ) 2 



6 



N. 



26 



Ss- 



2i) 



(40) 



32tt V v ) V $ 

We see from this estimate that the amplitude can grow quite steeply, especially for larger 
5, and the suppression due to the small coupling m v /v can be overcome by the growth in 
the number of states. Thus, for the theory to be unitary up to ^fs ~ M*, Eq. implies 
a lower bound on 1/R. 

We have performed a detailed calculation by taking into account the KK masses. We 
perform a change of variables to write n in terms of m^ n \ c.f. Eqs. (J22j) and (|23|) and use 
the continuum approximation to convert the double sum in Eq. (|3H|) into a double integral 
to get 



Im(^) 



-2 D 2S I dm (n) 



167T ^ . /0 



dm (m) [m (n) m (m) ] <5 1 Im(A^ (nm) ). (41) 
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FIG. 2: Unitarity bounds for 5 = 3. The theory is unitary if the solid line in (a) corresponding 
to the given 1/R does not cross the dotted line (Im(ao) = 1) at a yfs less than M*. The unitarity 
constraint on the 1/R vs. y/s plane is shown in (b). 



We use FeynCalc |22| to calculate M^ nm \ and we get 

M {nm) = K 2 hB {s,mW\mW 2 ) + 2B {t,ml,ml) +2B {u,ml,m 

+ f°(s) C (m 2 h , ml s, m^ 2 , m 2 , m'"' 2 ) 

/ 2 c (t) C (m 2 h , ml t, ml, m (m) \ m 2 v ) + (t -> u) 

~fc(t) C (ml ml t, ml, m^ 2 , m 2 u ) + (f - u) 
f\t) D (m 2 h , m 2 h , m 2 h , m 2 h , s, t, 2 , m 2 , m^ 2 , m 2 v ) + (t u 



(42) 



+ 
+ 

+ 



where B , C and D are the Passarino-Veltman (PV) scalar functions [23], and j{, f%, f£ 
and f d are defined as 

fi(s) = 4 (s + 2m 2 u + 2m {m) m u + 2m {n) m v + 2m (m) m (n) - 2m 2 h ) , (43) 
/ 2 c (t) = t + 2m (m) 2 + 4m„m (m) + 2m 2 - 2m 2 , 
/ 3 c (t) = f + 2m (n)2 + 4m u m (n) + 2m 2 - 2m 2 h , 

f d (t) = t ((m (m) - m (n) ) 2 - s) + 2 ((m (m) + m,) 2 - m 2 ) ((m (n) + m,) 2 - m 2 ) . 

The functions /f, /|, /| and / d are all real, and therefore, Im(.M ( - nm ' ) ) is given by the 
imaginary parts of the PV scalar functions. We present the calculation of the imaginary 
parts of the PV functions B , Co, D , in Appendix El To obtain Im(^4), we then perform 
a Monte-Carlo integration of Eq. (|41j) over m^ 1 ' and m^ m ' , and expand it into partial waves 
using Eq. (|3*5jl. We show the resulting Im(a ) in Fig. ^a) for 5 = 3, as a function of t/s for 



15 



different choices of 1/R (from 10 2 eV to 10 4 eV). For 5 = 3, the region in 1/R vs. that 
satisfies the unitarity bound, c.f. Eq. (jHEJ), is shown in Fig. Efb). We infer from these plots 
that for 5 = 3 and y/s = 1 TeV, we require 1/R> 250 eV, for unitarity not to be violated. 

Using the bound just obtained for 5 = 3, we can use Eq. (pfUj) to estimate the bounds for 
5 = 1,2. For 5 = 2 and y/s = 1 TeV, we find the unitarity bound, 1/R > 1CT 3 eV. For 
5 = 1, we find that the unitarity bound is very weak. 



V. CONCLUSIONS 

The Standard Model (SM) of high energy physics suffers from the gauge hierarchy prob- 
lem, which is the fine tuning required to maintain a low electroweak scale (Mew ~ 10 3 GeV) 
in the presence of another seemingly fundamental scale, the Planck scale (the scale of grav- 
ity, M p i ~ 10 19 GeV). Theories with n large extra dimensions have the potential to address 
the hierarchy problem as shown by Arkani-Hamed, Dimopoulos and Dvali (ADD) 0, 0| . 

Recent neutrino oscillation experiments H IE have suggested a small non-zero mass 
for the neutrino. It has been shown [ESE^I that if right-handed bulk neutrinos (SM gauge 
singlets) propagate in 5 of these large extra dimensions (5 < n), then the smallness of the 
neutrino mass can be explained. 

We considered such bulk neutrinos in the ADD framework where the fundamental scale 
in nature is postulated to be M* ~ 1 TeV. Theories with 5 > 2 naturally give the right order 
of neutrino masses with the coefficient of the Yukawa interaction for the neutrino being 
0(1). However, for 5 = 1, with the requirement that oscillation constraints be satisfied, this 
coefficient has to be O(10 -5 ), c.f. Table HJ making the theory unnatural. 

We performed the Kaluza-Klein (KK) expansion of the bulk neutrino field to obtain the 
equivalent four dimensional theory starting from a 4 + 5 dimensional theory. It is common in 
the literature to analyze theories with 5 = 1,2, while in this paper we considered 5 = 1,2,3. 
However, in the case of 5 > 2, we find, as expected, for some observables the sum over KK 
states can be divergent. We took an effective theory approach and cut the KK sum off at 
M*, where we expect new physics to regulate this divergence. 



We take the view, as others have 



121 ]. that the standard three-active-flavor analysis 
explains the neutrino oscillation data. (We do not address the LSND result.) This seems 
highly plausible given the Super-Kamiokande and the SNO neutral current data. Thus the 
oscillation into sterile KK neutrino states must be constrained to be smaller than the current 
precision in experimental data. 

Using the KK theory, we estimated the probability of the three active neutrino species to 
mix into sterile KK neutrino states. This probability is constrained by experiments such as 
CHOOZ and fits to the atmospheric oscillation data, leading to bounds on the radius (R) of 
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extra dimensions that right-handed neutrinos propagate in. We have compiled the bounds on 
1/R for 5 = 1, 2, 3 in Table. ITTl showing a strong bound for 5 = 3. For example, in the normal 
hierarchy of neutrino masses, we find 1/R has to be bigger than about 0.15, 1.5, 5600 eV for 
5 = 1,2,3 respectively. We also present bounds for the inverted and the degenerate mass 
schemes. The mixing probability is divergent when summed over the KK states for 5 > 2, 
with a mild logarithmic dependence on the physics that comes in at M* for 5 = 2, and a 
stronger linear dependence in the case of 5 = 3. We do not predict the bounds for 5 > 3, 
since the dependence on the high energy completion of the extra dimensional theory we are 
working with would be increasingly stronger for bigger 5. Our bounds have an additional 
source of uncertainty from using the continuum approximation rather than performing a 
discrete sum over the KK states. 

We also present bounds on 1/R coming from maintaining perturbative unitarity in the 
theory, in the Higgs-Higgs scattering channel. The bound is due to the growth in the 
large number of KK states that contribute in the intermediate state as the collision energy 
increases. The unitarity bound is also quite strong, especially for 5 = 3, as shown in Fig. |2l 
albeit not as strong as the bound derived from oscillation data. 

It is interesting to ask what signatures of bulk neutrinos there might be at a high energy 
collider. Since the only coupling to SM fields of the bulk neutrino is through the Yukawa 
coupling, one promising channel is the production of the KK states in the final state in 
association with the Higgs. Even though the coupling of the Higgs to the neutrino is much 
suppressed, the rate can be enhanced to measurable levels owing to the large number of KK 
states that can be produced in the final state. This is the subject of our forthcoming paper. 
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APPENDIX A: 2-COMPONENT DIAGONALIZATION 

We can arrive at the mass eigenvalues obtained in Sec. Ill Bl using two component Weyl 
fields and putting the neutrino mass matrix in a symmetric form. This is achieved by making 
the field definitions [8], 




(Al) 



We can then rewrite Eq. (0) as, 



£( 4) = _1J| ^ N (n) N (n) _ M (n) M (n) + ^ 



m 



v 



v 



h u R u L + Y j {N^u L + M^u L )+h.c.[A2) 



n 
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Defining v T = (u L u R . . . M^) and with (h) = v, we get the neutrino mass 

term, 



(A3) 



where the neutrino mass matrix M. is given by, 
/ m v m u m v 



M 



m,. 



\ 



m v 
m v 



l_ 

R 



m, 



N 
R 



V 



• (A4) 



R ' 



The alternating sign on the diagonal reflects the Dirac nature of masses. We have already 
noted that for 5 > 1, the state with mass ±^ at the n th level has a degeneracy d^. 



The characteristic equation, Det(.M — A 1) = is, 



N 



n (n 2 -(^) 2 ) 
\n\=i 



2\ d A 



N 



dt 



(XR) 2 - (m u R) 2 + 2(XR) 2 (m v R) 2 V - 

l*[=i 11 v 7 



0. (A5) 



where, similar to n, is shorthand for (ki, k$). In order to respect neutrino oscillation 
experimental bounds, we have to restrict ourselves to m v R = £ <C 1, and in the analysis 
that follows, we will keep terms only up to 0(£ 2 ). 

The lightest neutrino mass given by the lowest eigenvalue of Eq. (|A5[) . in the limit 



pa 



(A6) 



The heavier masses are as follows: 
For 5 = 1, we have d n = 1 for all n and this case was considered in Ref. J8j. For the KK 
level with d n = 1, in Eq. ()A5jl . Ai? = n is never a solution, and only the second of the two 
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factors in square brackets is relevant. We find the eigenvalue \ n R = ±n(l + £ 2 + £ 4 ) and 
the corresponding eigenvector, 

-L 1 , . ™ ™ „, s T 



(±A„-fc/fl) (±A n +fc/R) 



(A7) 



where the superscript T on the row vector represents the transpose operation, and the 
reciprocal of the normalization constants are, 



1 

No 



d;, 



— = £ 2 d~ 



However, at the n th KK level, if d n > 1 we can infer from Eq. (jA5|) that there are two classes 
of eigenstates: 

• one with eigenvalue ±^r [1 + (£ 2 + ^ 4 )d n } with eigenvector as shown in Eq. (|A7|) . 

• (d n — 1) states with eigenvalues and the corresponding eigenvectors having the 



first 2n entries zero, and with the structure, 



J n(2,...,d n ) 







v x 



V2d n 







(A8) 



The eigenvalues found here agree with Eqs. (^|) and (J2BI)- 

The matrix V that diagonalizes the mass matrix Eq. (|A4|) have the e n as columns. Neu- 
trino oscillation from an active state to a sterile state computed using this matrix V agrees 
with that found using Eq. (|2U|) . 



APPENDIX B: IMAGINARY PART OF PASSARINO-VELTMAN FUNCTIONS 

The unitarity of the S— matrix requires that twice of the imaginary part of any scattering 
amplitude can be written as 

2lmM(a ^b) = J2 [ dIL f M*{b -> f)M{a -> /), (Bl) 

in which the sum Y2f runs over all possible final-state /, and Uf represents the phase space 
of final-state /. We use the Cutkosky cutting rules to calculate the imaginary part of the 
Passarino-Veltman scalar functions B ,C , D which were encountered in Sec. IIVI 

The Feynman diagrams (Aii) and the corresponding cut diagrams which are re- 

lated with each other as in Eq. (|B1J) . are shown in the upper and lower parts of Fig. 01 
respectively. The label i = 2,3,4 represents the two-point, three-point and four-point dia- 
grams, respectively. In general the N-point one- loop scalar integral in D— dimension is 

T N (Pu • • • , Pn-i, rajv-i) = ^— ^ f d D g 1 , (B2) 

m J L> L> 1 ■ ■ ■ P N _i 
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with the denominator factors 

V = q 2 -ml, T>i = (q + p { ) 2 - m 2 , with i = 1, • • • , N 

originating from the propagators in the Feynman diagram. The ie part of the denominator 
factors is suppressed. Here we adapt the conventions of Ref. 



1. Two-point integrals, Bq function 

The scalar two-point diagram (A4 2 ) in Fig. Eta) is given by 

. j f d 4 q 1 1 , 2 2 2 , 

Ado = / ~, 7-7 — — = ^Bn[p ,777.1,777.9), 

J (2tt) a V V 1 16tt 2 ^ ' u 2h 
from which we obtain the imaginary part of _B , by applying the Cutkosky rules, as 

Im(S Q ) = 8vr 2 x 2Im(A^ 2 ) 
= 8tt 2 x M' 2 

7 2E 1 2E 2 V ; 
= tt/T, (B3) 

where (3* is defined as 

(3* = \\/(P 2 ~ ( m i + m 2) 2 ) ip 2 - (mi - m 2 ) 2 ). 
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2. Three-point integrals, Co function 

The scalar three-point diagram (M3) in Fig. OJb) is given by 

from which we obtain the imaginary part of Cq function as 

Im(Co) = 8tt 2 x M' 3 

1 _ f 1 



4 7 (£4 - pi) 2 - m|' 
In the cm. frame, with the scattering angles 8,<ft,6*,<j)* defined in Fig. 0J the 4— momenta 



is 



Pi = ^(l,/?sin0,O,/3cos0), 
p 2 = ^(1, -P sin 0,0,-/3 cos 9), 
fci = ^-(£; 1 ,^*sinrcos0*, ( 5*sinrsin0*,^*cosr), 
fc 2 = ^-(E 2 ,-/?*sinr cos0*,-/?*sinr sin0*,-/?*cosr), (B4) 
with s = p 2 = (pi + P2) 2 , pi = v\ — m a> ^1 = m i> ^2 = m 2 an d 



= 



_ ^/(g - (mi + m 2 ) 2 ) (s - (mi - m 2 "pj 

s 

- s + m 2 — m 2 

-c/l = , 

S 

E 2 = g + m2 " m \ (B5) 

s 

After a straightforward calculation, we find the imaginary part of Cq to be 

ird* 1 f + 33* 

where / is defined as 

s — m 2 — mj — 2m 2 + 2m\ 
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FIG. 4: Kinematics of the cut diagrams 



3. Four-point integrals, Dq function 



The scalar four-point diagram (A4a) in Fig. |3^c) is given by 
d A q 1 1 



M 4 



(2vr) 4 VqV{D 2 V z 16tt 2 
from which we obtain the imaginary part of D function as 



7-1/2222 j. 2 2 2 2\ 
D {P 1 ,P 2 , P31 S > *j m H m 3' m 2> m 4)' 



Im(£) ) 



Using the same kinematics defined in Eqs. (jB4|) and (jB5 



8tt 2 x Af 4 

If 1 1 

4 7 (h -pi) 2 - m\ (ki - p 3 ) 2 - mj' 



lm{D 



2ttP* 1 



In 



[BP - 2AQ){Px + Q)+2K- 2P^K(Ax 2 + Bx + C)' 



2irp* 



arcsm 



Px + Q 



2K + (BP - 2AQ)(Px + Q) 
P(Px + Q)VB 2 -AAC~ 



x=l 



2np* 



2P VAx 2 + Bx + C 
{BP-2AQ){Px + Q)_ 



x=l 



x=-l 



(for K > 0), 
(for K < 0), 
(for K = 0), 
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where the variables A, B, C, P, Q and K are defined to be 

A = (3 2 /3* 2 , (B8) 
B = -2/ p (3* cos 6, 
C = f - f3 2 (3* 2 sin 2 6, 

p = /3/r, 

Q = -f, 

K = AQ 2 - BPQ + CP 2 , 

and / is defined as 

s — m\ — mf, — 2m 2 + 2m\ 
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